Simulation of indivisible qubit channels in collision models 
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A sequence of controlled collisions between a quantum system and its environment (composed of 
a set of quantum objects) naturally simulates (with arbitrary precision) any Markovian quantum 
dynamics of the system under consideration. In this paper we propose and study the problem of sim- 
ulation of an arbitrary quantum channel via collision models. We show that a correlated environment 
is capable to simulate non-Markovian evolutions leading to any indivisible qubit channel. In par- 
ticular, we derive the corresponding master equation generating a continuous time non-Markovian 
dynamics implementing the universal NOT gate being an example of the most non-Markovian quan- 
tum channels. 

PACS numbers: 03.67.Lx,03.65.Ta 



I. INTRODUCTION 

It is one of the basic postulates of quantum physics that 
the time evolution of closed quantum systems is governed 
by the Schrodinger equation [l[ . As a result for any time 
interval of length r > a state of a quantum system de- 
scribed by a density operator g is transformed according 
to a unitary conjugation 



U tQ UI , 



(1.1) 



where U T = e~»i , H is the (for simplicity time- 
independent) Hamiltonian of the system and h is the 
Planck constant. Moreover, it is natural to assume that 
U T = UtU s , where r = t + s and t, s > are arbitrary du- 
rations of short time intervals that sum up to the whole 
time interval r. However, if the system is open (inter- 
acting with its environment) such simple division of time 
evolution into incremental time steps does not necessarily 
possess a clear meaning. It is a relatively recent discov- 
ery that certain non-unitary quantum evolutions are 
indivisible into shorter ones. In order to make this state- 
ment precise let us introduce a mathematical framework 
modelling open system dynamics over finite time inter- 
vals. 

Let us assume that initially the system and its envi- 
ronment are not correlated, so that that the combined 
system-environment initial state can be described as f><8>£ 
(such an assumption is quintessential in order to preserve 
linearity of the dynamics - see eg Ref. Q and refer- 
ences therein). Then the evolution of the joint system- 
environment state is governed by the Schrodinger equa- 
tion and the state transformation over the time interval 
[0, r] is described by a map 



q i-> £ T [g] = tv env [U T g <g> 0}] 



(1.2) 



where £ is the initial state of the environment and tr env 
denotes the partial trace over the environmental degrees 
of freedom and U T is a unitary transformation describ- 
ing both the system and the environment. Abstracted 



mathematical properties of £ T guaranteeing the exis- 
tence of such model for £ T are linearity (£(X + XY) = 
£(X) + X£(Y)) 7 complete positivity (positivity oi£ T ®X d 
for all d) and trace-preservation (tr[£ T [X]] = tr[X]). This 
result is known as the Stinespring representation theorem 
[i[ and the maps £ T with such properties we call quantum 
channels (see for example [f|). 

We say that a quantum channel £ is divisible if £ = 
£io£ 2 , where £\, £i are both non- unitary quantum chan- 
nels. If this is not the case, then we say the channel is 
indivisible. According to this definition unitary channels 
are indivisible. This may sound quite counterintuitive, 
especially, if we recall how unitary operations are decom- 
posed into gates in quantum computation Q . Indeed the 
gate decomposition is the crucial tool in the analysis of 
quantum complexity. But let us stress a difference be- 
tween the decomposition and the divisibility of channels. 
The primary task of the decomposition is to simplify 
the implementation of a multipartite quantum channel 
by means of elementary gates, thus, the decomposabil- 
ity captures the complexity of the process expressed in 
terms of number of uses of elementary gates. On the 
other side, the divisibility addresses more fundamental 
question whether a given process can be understood as 
a sequential concatenation of processes. If not, then it 
is natural to ask how to simulate such quantum chan- 
nels, in particular by means of some "continuous" time 
evolution. 

Dynamics of open quantum systems is modelled by so- 
called master equations that are usually derived under 
the assumption of Markovianity 0, [1] . Mathematically 
this means that one parametric set of channels £ t , being 
the solution of Markovian master equation (cquivalcntly 
Lindblad master equation), satisfies the semigroup prop- 
erty £ s o £ t = £ s+t for all s,t > and the initialization 
condition lim t _>o £t = Formally we can write £t = e ct , 
where C is known as the Lindbladian. We say a channel £ 
is Markovian if there exist the Lindbladian C and time t 
such that £ = e Cr . Otherwise, the channel is called non- 
Markovian. Let us note that all indivisible channels are 
non-Markovian. but there are divisible non-Markovian 
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channels @, Q . 

Simple collision models, in which the system's evolu- 
tion is modelled via a sequence of (weak) interactions of 
the system with sequentially selected particles from the 
environment, provide a natural playground for simula- 
tions of quantum dynamics. It is known (see for instance 
Refs. [HKIlj]) that these models can approximate (with 
arbitrary precision) any evolution governed by a Lindblad 
master equation [12l. Il3j|. i.e. implement any Markovian 
quantum channel. Our aim in this paper is to investi- 
gate whether a general quantum channel (in particular 
an indivisible one) can be simulated in the framework 
of simple collision models. Let us note that a related 
collision model was studied in Rcf. where the au- 

thors derived master equations for a family of correlated 
Markovian evolutions. 

This paper is structured as follows. In Section II, we 
define the framework of collision models and we intro- 
duce the concept of stroboscopic simulations. In Section 
III, we design a collision model for implementation of all 
indivisible qubit channels. Finally, we discuss our results 
and present conclusions of our investigation in Section 
IV. 



II. SIMPLE COLLISION MODEL 

By a collision model of a specific system dynamics 
we understand a sequence of interactions (collisions) be- 
tween the system and particles from the environment. 
Essential property of the model considered in this paper 
is that each particle of the environment interacts with the 
system at most once, while the environment particles do 
not interact between themselves. More formally, let Uj 
be a bi-partite unitary operator describing an individual 
collision between the system particle and the j-th particle 
of the environment. Collision model is a concatenation 
of unitary channels Ui0.03U n (see Fig.[T]). Let us denote 
by ui n the initial state of the environment composed of n 
particles and by g the initial state of the system. 
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FIG. 1: Schematic illustration of a simple collision model. 



Suppose g — > Qt = £t (g) is a solution of some Lindblad 
master equation. The question we will like to answer is 
"How to simulate the process £ t in the framework of a 
collision model?" Assume Us and £ determine the Stinc- 
spring dilation of £$. Then the sequence of concatena- 
tions of Us applied in each step jointly on the system 



and a "fresh" particle from the environment in the state 
£, results in a discrete evolution 

Q n =£%{q) 

= twKE/i • • • U n )(g ® e n )(Ui ■ ■ ■ C/ n ) f ] , (2.1) 

where Uj = Us <£> Ij and Ij denote the identity opera- 
tor on all particles of the environment but the jth one. 
As a result the described collision model stroboscopically 
simulates the continuous Markovian evolution, i.e. for 
all 5 > we have £ t ~ £ J s l with n = [t/S] ([x] denotes 
the nearest integer to x). The parameter 5 determines 
the quality of the (stroboscopic) approximation of the 
time-continuous (Markovian) dynamics £t for t > 0. 

It was shown in [l5| that if U is a controlled unitary 
operator, then the sequence of collisions simulates pure 
decoherence process, in which the diagonal elements of 
density operators (with respect to the so-called decoher- 
ence basis) are preserved while the off-diagonal ones van- 
ish exponentially. Setting U to be the partial SWAP in- 
teraction [l6|, [TtJ the system exponentially converges to 
the original state of the particles in the reservoir. The in- 
formation about the original state of the system is diluted 
into the correlations among environmental particles. Be- 
cause of these features the process is called quantum ho- 
mogenization and represents a quantum information ana- 
logue of thermalization process, in which the temperature 
is replaced by the concept of quantum state. 

In both examples described above it is assumed that 
the initial state of the reservoir is factorized, i.e. ui n = 
£® n . In this factorized setting, the channel £\ is called 
a generating channel of the collision process. Indeed, 
if this is the case, then the induced discrete dynamics 
£ n = £[ l fulfills the conditions of a discrete semigroup, 
i.e. £ n o £ m = £ n+m for all positive integers n, m. With- 
out loss of generality we may set £q = X. It is important 
to stress that such semigroup feature does not guarantee 
that the discrete dynamics stroboscopically approximates 
some Markovian continuous-time evolution. In other 
words, the generated channels £ n are not necessarily 
Markovian. For instance, by definition indivisible chan- 
nels ^indivisible can be elements of a discrete semigroup 
only as the generating channels, i.e. £ indivisible = £\- 

But what if a general initial state of the environment 
is allowed? Could we then find a way how to stroboscop- 
ically simulate (using the collision model) a continuous- 
time evolution leading to an arbitrary quantum channel? 

Let us formulate the problem: We say £ is a strobo- 
scopically simulated channel if for every n there exists a 
bipartite interaction U (between the system and a spe- 
cific particle from the environment ) and an initial state 
uj n of the environment such that £ n = £ , where 

£ n {g) = tawKPi • • • U n )(g ® w n )(Ui ■ ■ ■ , 

and Uj = U ® Ij. If the channel can be approximated in 
the introduced sense, then the collision models reveal the 
elementary features of Markovian dynamics even for non- 
Markovian channels. In particular, the process can be 
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realized in arbitrarily small (non-unitary) steps although 
the steps themselves are not described by valid channels. 
Our aim is to address the question which channels can 
be stroboscopically simulated? In what follows we will 
focus on implementation of indivisible quantum channels, 
which are, intuitively, the most non-Markovian examples 
of quantum evolutions. 



A. Evaluation of ||fj+i — £j\\ 

In this Section we will show that for the considered 
collision model the individual collisions induce arbitrarily 
small disturbances of the system. Assume n is fixed, i.e. 
77 = 7r/ (2n). After jth interaction the state of the system 
undergoes the transformation 



III. STROBOSCOPIC SIMULATION OF 
INDIVISIBLE QUBIT CHANNELS 

The existence of indivisible channels was observed for 
the first time in Ref. @, where the authors analyzed in 
detail qubit channels. In particular, they found that in- 
divisible qubit channels are unitarily equivalent to the 
subclass of Pauli channels 



£{q) = <lxO x Qo x + q y a y ga y + q z a z ga z 



(3.1) 



where q x ,q y ,q z are positive {q x q y q z ^ 0) and q x + q y + 
q z = 1. For instance, the choice q x = q y = q z = 1/3 de- 
fines the best quantum approximation f not of the uni- 
versal quantum NOT gate for a qubit J-not : Q i-> (I~o) 
[ill, [l9|]. The question is: how to simulate stroboscopi- 
cally some continuous time dynamics leading to these 
channels? 

Suppose environment is composed of three-dimensional 
quantum systems (qutrits). Then one of the possible 
dilations for indivisible qubit channels is the following: 
£{g) = tr cnv [C£ ® wC*], where C = a x <£> \x)(x\ + 
dy ® \y)(y\ + <J Z ® \ z )( z \ is a controlled unitary inter- 
action and lu satisfies the conditions (k\uj\k) = q k for 
k = x,y, z. We assume that |x), \y), \z) form an orthonor- 
mal basis of the qutrit's Hilbert space. Let us stress that 
C is not only unitary, but also Hermitian. Therefore 



ir)C 



cos 777 + 1 sin rjC defines a family of in- 



teractions. Moreover, each member of this family is a 
controlled-U operator, U v = ^2 k e lTI(Tk ® \k)(k\. 

Consider a collision model generated by collisions U n . 
The concatenation of n such collisions implements a 
global (system plus environment) unitary transformation 

Ur,,! ■ ■ ■ U v>n =Y d ^nr,a k ] ® \k® n )(k® n \ . (3.2) 



Suppose u>„ is the initial state of the environment com- 
posed of n qutrits such that (k® n \uj n \k® n ) = q k . After 
nth collision the system evolution reads 

£ n (e) = tr env [(C/ r(i i---C/ 7)j „)(e(8)a;„)([/^i---J7 ?)j „) t ] 
= X] qk ex P[ inr l^k} Q exp[-inr]a x ] , 



hence setting the strength of the interaction 77 = n/(2n) 
the resulting collision model simulates the implementa- 
tion of any indivisible qubit channel in n steps. 



k 

-sm( — ) > q k [a k ,g , 

rj — ' 

k 



2n 
qkVkgvk 



n 



(3-3) 



where [A, B] = AB — BA is the commutator of operators 
A, B. Let us denote by £ the target quantum channel, 
i.e. £ = £ n (we fix 77 = ir/(2n)) and define .F(-) = [F, ■] 
with F = i^2 k qk&k- Then the distance between two 
subsequent steps equals 



A 

where 



3+1 



£, 



\C J+1 (2-£)+D 1+1 F\\,(3A) 



Cj+i 



cos 

1 



(j + 1)tt 2 jir 



2n 



(J + 1)tt 

cos cos 

2 V n n 



2n 



U + 1)tt 



Dj+i = - j sin — 



n 



(3.5) 
(3.6) 



Let us note that the norm we have in mind here is the 
completely bounded norm (see for instance (20j). namely, 



sup \\(A ®I„)(A)||i, 

n,X:||X||i=l 



where 



tr • 



Using the identities cos 2 a = (l + cos2a)/2, cos a — 
cos/3 = -2sin(^±^)sin(^ 2 ^), and sina - sin/3 = 
2sin(^)cos(«±£) we get 

. (2j + 1)tt . tt 
Cj+i = — sin sin — , (3.7) 



D 



3+1 



(2 7 + 1)7T TT 

cos sin — 

n n 



(3.8) 



This allows us to conclude that the distance is bounded 
as follows: 



A < |C J+ i|-||I-f|| + |D i+1 |-||^|| 

\I-£\ 



TT 

< sin — 

n 



. (2j + 1)tt 
sin 



< K sin — ( 
n 



n 

for n 



cos ■ 



(2j + 1)tt 



00) , 



where A' < ||Z - f || + ||J"|| < 2 + ||J"|| < 00, be- 
cause for channels \\£\\ = 1 and \\F\\ < 00. In 



(3.9) 
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fact, as shown in Rcf. [20J for any linear qubit 
map H^ll < 2^sup|| X || 1=1 ||^(X)||i, hence, ||F|| < 

2V2su P||A . ||i=1 \\FX - XF\\! < 4V2^)ji < 4^||F||, 
where the relations ||FX||i < ||F|| • ||A||i and = 
sup^ HFVII/IMI < 2 were used. That is, K < 2 + 8V2, 
which allows us to conclude that the system's changes 
in individual steps can be made arbitrarily small as n 
goes to infinity. Therefore, in this limit the evolution is 
continuous. 



In the Bloch sphere parametrization the channels take 
the form of an affine 4x4 matrix. Define 3x3 matrices 



I = 



S 





(3.12) 



B. Master equation 

Replacing the integer parameter j in the expression 
for Sj by a continuous parameter t we formally define a 
one-parametric set of channels 

St — £ + cos 2 (at) (I — S) + sin(at) cos(at)F , 

where a = w/(2n) and linear maps S , F are defined in the 
previous paragraph. Let us note that £ t form the same 
set of channels irrelevant of the value of a. Moreover, 
this one-parametric set of channels is continuous in t and 
So = X. A collision model (determined by the value of 
a) stroboscopically simulates continuous-time quantum 
evolution given by £ t , where different values of a define 
the quality of the simulation. In this section we will 
derive the master equation generating £ t . 

Formally, the evolution of density operators is gener- 
ated by the first order differential equation 



dg t 
dt 



d£t 
dt 



t{0t) 



(3.10) 



where 



+[<TjX, a k }) 



(3-11) 



is the generator of the dynamics [L2J, [l3[ . Any such gen- 
erator defines a dynamics which is tracepreserving and 
if (time-dependent coefficients) hj,Cjk are real, then also 
the hermiticity of operators is preserved. If these parame- 
ters arc time-independent {C t = C) and matrix composed 
of entries Cjk is positive, then the generated dynamics is 
also completely positive and Markovian. In such case we 
can write £ t = e ct and C is the Lindbladian. 

For sake of simplicity let us illustrate the derivation of 
the driving master equation (£ t ) for the case of the target 
channel £ = £not- This channel transforms the Bloch 
vector r into f = —r/3, hence, implements a Bloch ball 
(shrinking) inversion. In this case the map T induces 
the Bloch vector transformation r — > u x r with u = 
(1/3, 1/3, 1/3). Let us stress that IF (I) = O, hence, the 
map is not trace-preserving and maps any operator into 
a tracclcss one. 



Then the one-parametric dynamics is given by the fol- 
lowing 4x4 matrix 



10 

x(t)I + a(t)A 



where x(t) = -|(4cos 2 at — 1) and a(t) = ^sin(2at). A 
direct calculation gives 



dS± 
dt 











1 o 

3 a \0 4sin(2at)7 + cos(2at)A 



(3.13) 



and 



£7 X 















-I 



3x(l) dct £± 



[x{t)A - a{t)S] 



(3.14) 



where det£ t = 3x(t)[a(t) 2 + x(t) 2 ]. Thus, for the gener- 
ator we get 



with 
b(t) 
c(t) 

d(t) 



dStp-x _ 
dt 1 ~ 















-a(t) 



9x(t) 



12a 



b(t)I + c(t)A + d(t)S 



1 



, (3.15) 



a{tf +x(t) 2 
a(3x(t)-l) a(t)[3x{t) + a(t)} 



(3.16) 



x(t)[ a (ty + x {ty 



a{t)[3a(t) - x(t)] 
9x(t)[a(t) 2 + x(t) 2 



(3.17) 
(3.18) 



Using the methods and formulas described in Ref. [ll| 
we obtain the non-Markovian master equation in the op- 
erator form 



dg t 
dt 



j 

-d(t)'Y]<TjQtak, (3.19) 



where H — a x + u y + a z . The time dynamics induced by 
the collision model is illustrated in Fig. [21 
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FIG. 2: The collision model simulating the continuous time 
evolution towards the universal NOT gate (shrinked Bloch 
sphere inversion). In particular, the transformation of the 
Bloch sphere (lines capture the time evolution of eigenstates 
of a z operator) is depicted for the time interval t £ [0,n]. 
For t = |n the channel £t in not invertible (det £t=2«/3 = 0) 
and at this time the Bloch sphere is mapped onto a two- 
dimensional disk. Let us note that images of eigenstates of 
CT Z operator are internal points of the disk. In fact, the whole 
disk is the image of pure states only. 



IV. CONCLUSIONS 

In this paper we open the question of continuous time 
(stroboscopic) simulation of quantum channels within 
simple collision models, i.e. via a sequence of interactions 
of the system with particles forming the environment. 
Using the environment composed of three-dimensional 
quantum particles we design a collision model simulating 
arbitrary indivisible qubit channel. Indivisible channels 
could be coined as the most non-Markovian ones if one 
quantifies the Markovianity of a channel £ as the maxi- 
mal number (n) of non-unitary channels £\, 0.03, £ n such 
that £ — £\ ■ ■ ■ £ n . The smaller the number, the more 
non-Markovian is the channel. Let us stress that there 
are several recent proposals for measures of Markovian- 
ity of continuous time evolutions [2l| - f23j . These measures 
can be applied not only to the derived master equation 
for implementation of £not, but can be also modified to 
the settings of discrete time evolutions, which are natu- 
rally generated by simple collision models. However, a 
more detailed analysis along these lines goes beyond the 
scope of this paper. 

Another important point we want to stress is the ne- 
cessity of initial correlations between the particles of the 
reservoir, which are introducing the memory mechanism 



present in any non-Markovian evolution. It is of inter- 
est to understand whether there is some deeper relation 
between the correlation structure of w„ and memory fea- 
tures of the system's dynamics (Markovianity). In or- 
der to answer this question one needs to understand the 
ambiguity of the stroboscopic simulation. For example, 
each convex decomposition into unitary channels induces 
a different collision model (as described in the previous 
Section) for the same channel £ . However, for this class 
of collision models there is no qualitative difference nei- 
ther in the initial states, nor in the derived evolutions. In 
the considered collision model, the correlations are rela- 
tively strong and all particles forming the environment 
are mutually pairwise correlated, but it could happen 
that there are qualitatively different collision models for 
which the structure of correlations is completely differ- 
ent, especially, it could be that particles entering the j'th 
and the fcth collision are initially uncorrclated if \j — k\ 
is sufficiently large. 

In summary, we have introduced the problem of stro- 
boscopic simulations of a general quantum channels. In 
this paper we reported the case study of qubit indivisible 
channels and we designed a collision model for imple- 
mentation of any channel from the family of random uni- 
tary channels (seeEJ. Since indivisible channels are the 
most non-Markovian, it is natural to conjecture that any 
channel can be simulated in the stroboscopic manner. In 
fact, it was shown in Q that nonunital qubit channels 
are infinitesimal divisible and can be approximated by a 
concatenation of Markovian channels that can be strobo- 
scopically simulated by factorized states of the environ- 
ment. Moreover, since unital qubit channels are necessar- 
ily random unitary and for them we have an explicit col- 
lision model (irrelevant of their Markovianity), we may 
conclude that collision models can stroboscopically ap- 
proximate any qubit channel. We believe the considered 
collision model deserves further investigation that finally 
results in a better understanding of non-Markovian fea- 
tures of general continuous time quantum evolutions and 
implementations of non-Markovian quantum channels. 
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Appendix A: Simulation of all random unitary 
channels 

In this Appendix we will extend the presented stro- 
boscopic simulation of qubit indivisible channels to any 
random unitary channel acting on the system of arbi- 
trary dimension. A channel is called random unitary if 
£ = V \ .iA ', . qj > 0, = 1 an d Vj are uni- 

tary operators. Suppose 1 < j < d and define a collision 
model generated by the interaction 

d 

£/ = £li>0W- 1/n , 

where V^ n are unitary operators such that (V^ n ) n = 
Vj. Assuming that initially uj = qj\j® n )(j® n \, we 



find (analogously as for the qubit case) that after kth 
collision 

£ k (g) = J2q J (V 1 k/n )g(V 1 k/ y, (Al) 
j'=i 

thus, £ n = £. Such construction works for any value of 
n and therefore we can conclude that arbitrary random 
unitary channel can be stroboscopically simulated. Let 
us note that this collision model defines a non-Markovian 
evolution also for Markovian channels. 
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